import numpy as np

A = np.mat([
    [23076.9, 0.000976562, -7692.31, -0.000152588, 1923.08, -0.000274658, -7692.31, 0, 1923.08, -6.10352e-05, -2884.62, -2403.85, -2884.62, 2403.85, - 2884.62, -2403.85, -2884.62, 2403.85],
    [0.000976562, 23076.9, -0.00012207, 1923.08, -0.000244141, -7692.31, -0.000152588, 1923.08, - 6.10352e-05, -7692.31, -2403.85, -2884.62, 2403.85, - 2884.62, -2403.85, -2884.62, 2403.85, -2884.62],
    [-7692.31, -0.00012207, 11538.5, 0.000488281, -2884.62, -2403.85, 0,
        0, - 2884.62, 2403.85, 961.539, 480.769, 961.539, -480.77, 0, 0, 0, 0],
    [-0.000152588, 1923.08, 0.000488281, 11538.5, -2403.85, - 2884.62, 0, 0,
        2403.85, -2884.62, -480.77, - 3846.15, 480.77, -3846.15, 0, 0, 0, 0],
    [1923.08, -0.000244141, -2884.62, -2403.85, 11538.5, 0.000488281, -
        2884.62, 2403.85, 0, 0, 0, 0, -3846.15, 480.77, -3846.15, -480.77, 0, 0],
    [-0.000274658, -7692.31, -2403.85, -2884.62, 0.000488281, 11538.5,
        2403.85, - 2884.62, 0, 0, 0, 0, -480.77, 961.539, 480.769, 961.539, 0, 0],
    [-7692.31, -0.000152588, 0, 0, -2884.62, 2403.85, 11538.5, 0, -
        2884.62, -2403.85, 0, 0, 0, 0, 961.539, 480.77, 961.539, -480.769],
    [-6.10352e-05, 1923.08, 0, 0, 2403.85, -2884.62, 0, 11538.5, -
        2403.85, -2884.62, 0, 0, 0, 0, -480.77, -3846.15, 480.77, -3846.15],
    # [1923.08, -6.10352e-05, -2884.62, 2403.85, 0, 0, -2884.62, - 2403.85, 11538.5, 0, -3846.15, -480.77, 0, 0, 0, 0, -3846.15, 480.77],
    [0, 0, 0, 0, 0, 0, 0, 0, 11538.5, 0,0, 0, 0, 0, 0, 0, 0, 0],
    [0, 0, 0, 0, 0, 0, 0, 0, 0, 11538.5, 0, 0, 0, 0, 0, 0, 0, 0],
    [-2884.62, -2403.85, 961.539, -480.77, 0, 0, 0, 0, -
        3846.15, 480.77, 5769.23, 2403.85, 0, 0, 0, 0, 0, 0],
    [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5769.23, 0, 0, 0, 0, 0, 0],
    [-2884.62, 2403.85, 961.539, 480.77, -3846.15, - 480.77,
        0, 0, 0, 0, 0, 0, 5769.23, -2403.85, 0, 0, 0, 0],
    [2403.85, -2884.62, -480.77, -3846.15, 480.77, 961.539,
        0, 0, 0, 0, 0, 0, - 2403.85, 5769.23, 0, 0, 0, 0],
    [-2884.62, -2403.85, 0, 0, -3846.15, 480.769, 961.539, -
        480.769, 0, 0, 0, 0, 0, 0, 5769.23, 2403.85, 0, 0],
    [-2403.85, -2884.62, 0, 0, -480.77, 961.539, 480.77, -
        3846.15, 0, 0, 0, 0, 0, 0, 2403.85, 5769.23, 0, 0],
    [-2884.62, 2403.85, 0, 0, 0, 0, 961.539, 480.77, -
        3846.15, -480.77, 0, 0, 0, 0, 0, 0, 5769.23, -2403.85],
    [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5769.23]
])

b = np.array([0, 0, 0, 0, 0, 0, 0, 0, 0, 1153.85,
              0, 576.923, 0, 0, 0, 0, 0, 576.923])

x= np.linalg.solve(A, b)

print ("Solution", x)
